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1. Introduction 

Entanglement is a property which is central to quantum mechanics, with bipartite 
entanglement being readily producible experimentally. There has been some progress in 
the production of tripartite entangled beams, with the entanglement often obtained by 
mixing squeezed vacua with linear optical elements PlEj- Other methods which create 
the entanglement using an actual nonlinear interaction have been proposed, using both 
cascaded and concurrent x*^^^ processes jHl IH El E] • In this article we investigate the 
fundamental limits to the achievable tripartite entanglement available from both the 
manipulation of squeezed states with beamsplitters and from a process which utilises 
concurrent nonlinearities. 

A two-mode system is considered to be bipartite entangled if the system density 
matrix cannot be expressed as a product of density matrices for each of the two 
modes. The definition of tripartite entanglement for three-mode systems is a little 
more subtle, with different classes of entanglement having been defined, depending 
on how the system density matrix may be partitioned [7j. The classifications range 
from fully inseparable, which means that the density matrix is not separable for any 
grouping of the modes, to fully separable, where the three modes are not entangled 
in any way. For the fully inseparable case, van Loock and Furusawa [Hj, who call this 
genuine tripartite entanglement, have derived inequalities which are easily applicable to 
continuous variable processes. In this work we will analyse two different Hamiltonian 
processes in terms of these inequalities, as well as in terms of three-mode Einstein- 
Podolsky-Rosen (EPR) P type criteria which we shall define. We note here that 
both these methods of detecting entanglement provide sufficient, but not necessary, 
conditions, so that one or the other may be more sensitive and useful in a given situation. 

2. Criteria for tripartite entanglement 

We will first describe inequalities which, if violated, demonstrate that a system exhibits 
true continuous variable tripartite entanglement. For three modes described by the 
annihilation operators dj, where j = 1,2,3, we define quadrature operators for each 
mode as 



sufficient to demonstrate continuous variable bipartite entanglement were developed by 
Duan et al ^2] and Simon jTTj. A set of conditions which are sufficient to demonstrate 
tripartite entanglement for any quantum state have been derived by van Loock and 
Furusawa [Hj. Using our quadrature definitions, the van Loock- Furusawa conditions 
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give a set of inequalities, 
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V^23 
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V{X^ - X2) + ViY^ + Y2 + ^3) >4, 
V{X^ - X3) + V{Yi + % + %) >4, 
V{X2 - X3) + V{Y, + + %) >4, 



(2) 



where V{A) = (A^) — (A)"^. As shown in reference jH], the violation of the first condition 
still leaves the possibility that mode 3 could be separated from modes 1 and 2, but this 
possibility is negated by violation of the second inequality. Starting with any one of 
the conditions thus shows that, if any two of these inequalities are violated, the system 
is fully inseparable and genuine tripartite entanglement is guaranteed. We note that 
genuine tripartite entanglement may still be possible when none of these inequalities is 
violated, due to the criteria being sufficient but not necessary. We also note here that 
the original van Loock-Furusawa correlations were written in a more complicated and 
general form, but for the symmetric systems we evaluate in this work, the form we have 
given is sufficient. 

3. Entanglement and Einstein-Podolsky-Rosen correlations 

The EPR argument was introduced in 1935 in an attempt to show that quantum 
mechanics not be both complete and consistent with local realism j^]. Schrodinger 
replied that same year by introducing the concept of entangled states which were 
not compatible with classical notions such as local realism ^21- In 1989 Reid |13j . 
and Reid and Drummond [H] proposed a physical test of the EPR paradox using 
optical quadrature amplitudes, which are mathematically identical to the position 
and momentum originally considered by EPR. Reid later expanded on this work, 
demonstrating that the satisfaction of the 1989 two-mode EPR criterion always implies 
bipartite quantum entanglement ^3]. It was also shown by Tan ^H] that the existence 
of two orthogonal quadratures, the product of whose variances violates the limits set by 
the Heisenberg uncertainty principle (HUP), provides evidence of entanglement. Tan 
demonstrated this in the context of teleportation, with the outputs from a nondegenerate 
optical parameteric amplifier (OPA) mixed on a beamsplitter. In this article we extend 
Reid's original approach, based on an inferred HUP between two quadratures, to the 
case of tripartite correlations, where quadratures of three different optical modes are 
involved. Just as the definition of tripartite entanglement is more complex than that 
of bipartite entanglement, with different classes of entanglement having been defined 
depending on possible partitions of the system density matrix we find that there 
is more than one way to define EPR correlations for a system exhibiting tripartite 
entanglement. We have used two of the possible methods in a previous publication [Sj, 
while Bowen et al have defined a third f7j. In this work we will give formal proofs 
that the correlations we defined previously also serve to demonstrate the presence of 
tripartite entanglement. 
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4. Tripartite entanglement from EPR correlations 

Genuine tripartite entanglement is verified if we rule out any bipartition of the density 
matrix (p), which is to say that the full system density matrix cannot be expressed in 
any of the following forms 

P = Y.Pr''p'^^ P = Y.PrPr''^ P = Pr"" Pr ■ (3) 

r r r 

If these factorisations are ruled out then so is the fully separable form p = J2r Pr^Pr Pr ■ 
There are two forms of the criteria that we need to consider, arising from one and two 
mode inference. 



4-.1. Experimental two mode inference 

The first inference scheme that we will use to prove tripartite entanglement and EPR 
correlations involves using experimental observations of two modes to infer properties 
of a third. The proof for this scheme is essentially the same as for the original bipartite 
entanglement result of ^^1; nevertheless, the proof exposes a freedom in the derivation 
that we wish to draw attention to, so we reproduce it here. We consider the separation 
of the density matrix in the form 

P = j:PrPr''- (4) 

r 

We now introduce the operators and y", {a G {A,B,C}) with = 2i. We 

see that the conditional probability of result x^ for a measurement of x"^ at A given 
a simultaneous measurement of x^ at B and cf" at C with results xf and xf is then 
P{x^\xf ,xf) = P(x^, xf , xp)/P(xf , xf ), where, assuming separability, 

P(x^,xf ,xf) = J2PrPrix^)Prixf,xf). (5) 

r 

Here Pr(x^) = (x'^lp^^lx^), and Pr(xf,xf) = (xf , xf Ipf*" |xf , xf ), where |x") are 
eigenstates of x°. Furthermore, normalisation of the density matrix requires YlrPr = 1- 
We may now use the measurements of x^, x^ to infer, with some uncertainty, the 
value of x^. The mean of the conditional distribution is then 

Mf^ ^ E P(.1.f , = Y. ""pj^lf^ {-\, (6) 

where {x^)r = ^^.a Pr(x'^)x'^. 

The variance Afx^ of the distribution P(x"^|xf,xp) is 



Now the mean-square "^^a Prix ){x — d) is minimised by the choice d = (x )r, so 
that 



Tripartite entanglement and three-mode EPR correlations 



5 



where (T^(x"^) is the variance of over the distribution Pr{x^). Since the estimate for 
the result x^ may not be optimal, we also define an error for the estimate A^^jg^jX"^, so 
that, after averaging over the results xf and xp, we find 

^lf,esti^ > ^IfX^ > E ^(^^- ^«'')^'^^ = T.Pr^rix% (9) 

and similarly A'^^ji)^ > '^^Pr^'riv'^)- Combining the two results and using the Cauchy- 
Schwartz inequality gives 

AL,xML/y^ > ypra'A^"") J2p,^',iy^) (10) 



> 



J2pr(TAx^)aAy' 



:iii 



For any p^, the uncertainty relation takes the form cr^(x^)cr^(?/^) > 1, so that the 
assumed bipartition of the density matrix (jlj) implies 

^IfS^^'^iy > 1- (12) 

The experimental observation of the EPR critierion A'^^j-x^Aj^j-y^ < 1 rules out the 
bipartition Tripartite entanglement is verified by ruling out all such bipartitions. 
The simultaneous experimental observation of the three criteria 

AlfX^Al^y^ < 1, (13) 
AlfX^Al^y^ < 1, (14) 
Alfx'^Alfy'^ < 1, (15) 

is sufficient to confirm tripartite entanglement. Note that the exact form of the 
expression involving xf , xf used in the inference (in this paper we deal with expressions 
of the form xf ± xc) does not enter into the derivation. We see that in fact this 
separability measure is entirely independent of the way information about the remaining 
subsystem is used to infer properties of a single mode. To be specific, if an A^-mode 
system density matrix is separable in the form p = J2r P^pl''"''"^'^'^^"'"^ , then, regardless 
of the way information from handled the inferred variances for mode k 

will satisfy the uncertainty relation Af^^x^Af^^y^ > 1, and the EPR criteria for this kind 
of separability follows. However, when > 3 there are additional forms of separability 
to be ruled out if genuine A^-partite entanglement is to be confirmed. In this work we 
will focus solely on tripartite entanglement. 

4-2. One mode inference 

The alternative scheme uses information about one mode to infer the combined 
properties of the other two. We consider the expression of the density matrix in the 
form 

P = J2p^''p?- (16) 
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The conditional probability of results x"^ and for measurements of and x^ at A 
and B given a simultaneous measurement of x^ at C with result is P( 
P{x^, x^ , xf) / P{xf), were, assuming separability 

P(x^,x^,xf) = ^p.P.(x^,x^)P.(xf). (17) 

r 

Here Pr{x^,x^) = (x^, x-^lp^^jx^, x^), and Pr(xf) = (xf |p.^|xf). 

The measurements of x*^ are used to infer, with some uncertainty, the combined 
quadrature operators x^^ = x"^ ± x^, y^^ = ± y^. The mean of the conditional 
distribution is 

l^f= E Pix^x^\xf)ix^±x^) = J2P-^^01{x^±x%., (18) 



where (x^ ± x^)r = ^^.a ^.b Pr{,x^, x^)(x^ ± x^). 



The variance Afx^^ of the distribution P(x'^, x-^|xp) is 

A.M" = E^^#r 5: P.(a:^x-)(x^±x--^:^^)l (19) 

Now the mean-square '^^a Pr{,x^ ± x'^ — c/)^ is minimised by the choice 
d = (x^ ± x^)r, so that 

> J^^^!^:^^^ Pr{x\x^){x^±x^-{x^±x%f (20) 

r- ^^^^ > 

= Jl'-rffr^ri^'^ ^21) 

where a'^{x^^) is the variance of x^ ± x'^ over the distribution Pr(x^,x^). Since the 
estimate for the results x^^ may not be optimal, we also define an error for the estimate 
^L/.est^i'^' so that, after averaging over the results xp, we find 

^lf,esAr > AL/^r > E ^(^f ) A?^r = (22) 

and similarly Aj^j-y^^ > YlirPr(^rhj±^)- Combining the two results, and using the 
Cauchy-Schwartz inequality gives 

Ai,xrAi,yr > Y.pr<yi{^i") T^p^M") (23) 



> 



(24) 



We now note that for any , the uncertainty relation for the combined quadrature 
takes the form a^{x^^)a'^{y^^) > 4, so that the assumed bipartition of the density 
matrix ()l(j|l implies 

AL/^r AL/^r > 4. (25) 
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The experimental observation of the EPR criterion A^^^x^^A^^^^^^ < 4, as given in 
(jnH|) . imphes inseparabihty. Genuine tripartite entanglement is verified by ruling out all 
such bipartitions. The simultaneous experimental observation of the three criteria 

AL;X^^AL,^r < 4, (26) 
^l,xi^^l,yf < 4, (27) 

AL/:^r AL/^r < 4, (28) 
is then also sufficient to establish tripartite entanglement. 

4-3. Practical criteria 

In practice one usually has access to certain moments of quadrature variables, in 
particular, the elements of the covariance matrix. It has been shown that optimised 
linear inference based on the covariance matrix is never better than knowing the 
full conditional probablity distribution for modes of interest. Consequently, if EPR 
inequalities are violated using linear inference then the more exact criteria developed 
above certainly are. Hence we can find sufficent conditions which may be of great 
practical value, in that experiments have been designed to measure them 



4-3.1. Two mode inference In this case we make a linear estimate Xi^est of the 
quadrature Xi for the mode i from the properties of the combined mode j + k, so 
that, for example 

Xi^est = a{Xj ± Xk) + c, (29) 

where a and c are parameters which can be optimised, both experimentally and 
theoretically fSl^l- It has been shown ^3] that this corresponds to minimising the 
variance 

KT/(X, - a{X, ±Xk)) = {{Xi - a{X, ± Xk)?) - {X, - a{X, ± Xk))' (30) 
with respect to a. The minimum is achieved when 

_ v{xa,±x,) 

vix,±x,) ■ ^''^ 

In the above V{A,B) = {AB) — {A){B). Defining the optimal inferred variance for X^ 
as V'^i{X;) = V::/{X,)\a=a^.^^, we obtain 

V{X,,X,±Xk)Y 

r"/(Xi) = V(XA - ^ . . (32) 

^ ' ^ ' V{X,±X,) ^ ' 

We follow the same procedure for the Y quadratures to give expressions which may be 
obtained by swapping each X for a y in the above to give the optimal inferred estimate 

V^^"/(r,) = V(Yi) - ^ ^ ^ — (33) 
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A demonstration of the EPR paradox can be claimed whenever it is observed that 
V^g/ {Xi)V^^^ {Yi) < 1, and we have shown that such an experimental outcome is possible 
whenever the theory predicts 

V'''^{Xi)V'''^{Yi) < 1. (34) 

Following on from what was proven above, this demonstration for the 3 possible values 
of i is then sufficient to establish tripartite entanglement. 

4.4- One mode inference 

In this case one measures Vinf{Xj ±Xk) = V^g/{Xj ±Xk ^ (^^i)\a=an^in- Linear inference 
leads to the expression for the optimal (minimum) variance in the inferred quadrature 
Xj + Xk 

V{Xi,X,)±V{X,,Xk)Y 
r"^(X, ± X,) = ViX, ± X,) - ^ ^— ^, (35) 



V{Xi 

which is merely a different form of the expression given in Ref . j3] , 

v{x,,x,±Xk)y 

V'^f{X, ± Xk) = V{X, ± Xu) - ^ ^. (36) 



V{X, 

The same expressions hold for the Y quadratures, and it is then straightforwardly shown 
that the HUP requires that 

V{X,±Xk)V{Y^±Yu)>^. (37) 

There is a demonstration of this three mode form of the EPR paradox whenever 
V^^^ {Xj ± XkjVlst {Yj ± Yfc) < 4, which is predicted to be possible when 

V'^^^i^Xj ± Xk)V'''^{Yj ± Yk) < 4. (38) 

As above, this demonstration for the 3 possible combinations also serves to establish 
inseparability of the density matrix. 



5. Entanglement via beamsplitters 



It is simple to show that one quadrature squeezed state, with squeezing parameter r, 
mixed on a beamsplitter with a vacuum input, results in a bipartite entangled state with 
a value of 2(1 + e"'') for the Duan criterion ^0], where a value of less than 4 represents 
bipartite entanglement. If an amplitude squeezed state is mixed with a phase-squeezed 
state on a beamsplitter, both with squeezing parameter r, this gives a value of 4e~^ for 
the same criterion. In this section we will quantify one possible way in which tripartite 
entanglement may be obtained, using squeezed states obtained from individual x*'^^ 
processes, which are subsequently combined utilising beamsplitters. 

A schematic of an apparatus which has been used by Aoki et al j2] to produce 
tripartite entanglement by this method is given in figure showing the three optical 
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0P02 




0P03 



Figure 1. The scheme which mixes three squeezed states on two beamsplitters. 



parametric oscillators (OPO) and the two beamsplitters used. They experimentally 
measured continuous variable tripartite entanglement, obtaining values for the criteria 
of (0) which were just above 3 using our quadrature definitions. In what follows we will 
first assume that OPOl produces an ideal minimum uncertainty squeezed state of the 
Y quadrature and 0P02 and 0P03 produce minimum uncertainty states squeezed in 
their X quadratures, all with squeezing parameter r. With the beamsplitter BSl having 
reflectivity yU and BS2 having reflectivity we may write expressions for the output 
operators hj in terms of the input operators dj as 



which allows us to calculate all the required output correlations in terms of the variances 
of the OPO outputs. The required variances are 



V{X,,) = {l-^)V{X,,)+^^V{X,,), 

V{X,,) = /i(l - u)V{Xa,) + (1 - /i)(l - l^)V{Xa,) + l^ViXa,), 
V{X,,) = fluViXaJ + iy{l - fl)V{Xa,) + (1 - ^)V{Xa,), 

V{Y,,) = {l-fi)V{Y,J + fiV 

= - lyWiY^J + (1 - /x)(l - u)V{Y,,) + uViYj, 
V{%,) = iiuV{YaJ + y{l - li)V{Ya,) + (1 - u)V{Ya.,), 



V^l - Ai fll + v^Ti ^2, 



Vyu(l - 1^) a-i - a/(1 - /u)(l - v) a2 + as. 



(39) 
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12 J 



vix,,,x,,) = - [fiViXaJ + (1 - ti)v{Xa,) - v{X 



(13/ 



(41) 



'u 1 - v) 



It is straightforward to see that the modes represented by hi and can be entangled, 
with the Duan criterion giving 



V{Xi,, - + V{%, + no) = 4 coshr - sinh 



(42) 



For the case of /i = 1/2 this simplifies to 4e~'', which is the well-known result for two 
ideal squeezed states mixed on a 50 : 50 beamsplitter. In the unbalanced case we will 
consider here, however, ^ = 2/3, and we find 



4 cosh r — sinh r, 



(43) 



(44) 



ViX,,-X,,) + Vin,+Y,,: 

which exhibits a maximal violation for 

1 / l + 2v/2/3 ^ 
r = - log — — ^ 1.76. 

2 ^\^l-2v^/3y 

It is of interest to note that the bipartite entanglement between these modes then 
disappears for r a little greater than 3, while, as we will show below, the violation of 
the tripartite entanglement inequalities continues to increase with r. 

Our idealised model shows that the violation of the van Loock-Furusawa inequalities 
increases with the degree of squeezing, so that, for minimum uncertainty squeezed 

St cltjGS, 



V{Xh^ — Xb^] 
V{Xh^ - Xb.^] 

V{Xb,-Xb,) 



l + /ii/ + 2V/i(l-/i)(l -u) 



1 + ji — liv + 2\/ /iz/(l — /i) 

2 - /i + 2/iA/z/(l - V 



2v//u(l-^)(l-z/) 
1 — fl + flU — 1\/ iiv{\ — /i) 



e', 



(45) 



V(n, + n, + ^3)= 1 + 2 W^v{\ - z/) + v//i(i - /t^)(i -y)^ V/i^^(l - /i) 



/i 



v/z/(l - z^) + V/^(l - /i)(l - /^) + V/iz^(l - /i)) 



e . 



With /i = 2/3 and z/ = 1/2 as in Aoki et al |2j, we find for the van Loock-Furusawa 
criteria of (j2I) 

5e-", (46) 



v^i2 = y^ 
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CO (units of y) 



Figure 2. The van Loock-Furusawa spectral correlations corresponding to (j2Jl for the 
Aoki scheme, with /i = 2/3, v = 1/2, and outputs calculated via a linearised fluctuation 
analysis of the standard OPO equations, with 7a = 7b = 1 and k = 10~^. The three 
curves are for different ratios of the pumping rates to the critical threshold pumping 
rate. All three correlations are equal for these parameters. The value Sij = 4 defines 
the upper boundary for true tripartite entanglement. 



indicating clear violation of Q and therefore genuine tripartite entanglement for 
r > 0.23. Note that, with the quadrature definitions we use, the reported results of 
Ref. j2j are all approximately 3.2. 

We may also now calculate the three-mode EPR correlations of section El for this 
idealised case. Again with /i = 2/3 and v = 1/2, we find for the inference of a combined 
quadrature via measurements on a single quadrature 

V^-/(n^+nj =;z7^, (47) 

so that 

which is less than 4 for any finite value of squeezing and hence exhibits three-mode EPR 
correlations and entanglement. The van Loock-Furusawa inequality is not violated for 
r < — log 0.8 ~ 0.22, which is not in contradiction with the EPR criteria as both provide 
sufficient, but not necessary, conditions. In this particular case the EPR correlation is 
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Figure 3. The van Loock-Furusawa correlations as for figure El as a function of tiie 
ratio of the pumping rates to the critical threshold pumping rate. The solid line is the 
smallest value at any frequency, while the dash-dotted line is the correlation at zero 
frequency. All three correlations are equal for these parameters. Note that values in 
the immediate vicinity of e/ec = 1 are of limited value in the linearised analysis used 
here. 



more sensitive to the presence of entanglement. For the inference of a single quadrature 
from the sum of the other two, we find 

= ^ , V'^'HY,) = , (49) 

^ ^ 2e^ + e-'- ^ ^ 2e-'- + e'- ^ ^ 

which gives the product 

v-Hxi^v-'iK) = j^^- (50) 

It is readily seen that this product falls below 1 for any finite r and hence gives another 
demonstration of EPR correlations for this system. 

In a more realistic analysis we must consider that OPOs do not produce minimum 
uncertainty squeezed states with a monotonically increasing squeezing parameter ^\ 
12011^, but exhibit quite different behaviours above and below the oscillation threshold. 
Here we consider output quadrature amplitudes at a particular frequency shift u, so 
following Collett and Gardiner we define the associated spectral quadratures Xj{uj) and 
Yj{uj), j = 1,2,3, and the associated spectral variances. We will consider an OPO where 
d and b represent the signal and pump modes, with respective cavity dampings 7a and 
76. With the classical pump represented by e and the effective nonlinearity by k, we 
find that there is a critical threshold pump value ec = '^alb/i^, below which the signal 
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Figure 4. The EPR spectral correlations, 5'™-'' {Xi)S'^'^^ (F^), for different ratios of the 
pumping rates to the critical threshold pumping rate. All parameters arc the same as 
in figure 121 



mode is not macroscopically occupied. The output spectra for this OPO system are well 
known ^Hj, with the below threshold spectral variances for mode a being 



1 + 



(51) 



while above threshold they are 



(27,76 -2Ke + a;2)2 + 72^2' 



47a (7fe 



(52) 



(^2_2/s:e)2 + (27a + 7,)2a;2- 
The output field d exhibits squeezing and three such outputs of three OPOs are used as 
the inputs di, 0,2, 0.3 of figure [H We note here that, as the results (jSH), are derived 
using a linearised fluctuation analysis, they are not valid in the immediate region of the 
threshold. Given this caveat, in figure El and figure El we display the results of using 
these spectral variances in the expressions for the van Loock-Furusawa correlations. In 
these figures, the Sij{uj) are the measurable output correlations which correspond to the 
Vij of 0, so that Sij{uj) < 4 implies genuine tripartite entanglement. It can readily 
be seen that the potential violation of the inequalities available from this system is 
much stronger than that measured so far by Aoki et al j2] and also that there is a large 
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violation far above threshold. Following Reid and Drummond it is also possible 
to use these expressions as inputs to calculate the EPR correlations of section El in the 
spectral domain. In figure ID we show results for the inference of one quadrature from 
a combination of the other two. In this value of S'''f{Xi)S'''f{Yj) < 1 indicates 

EPR correlations, and genuine tripartite entanglement. The results for inferring the 
combined quadratures from the single ones are found in this case by multiplying these 
results by 4, replacing, for example the with 5™-^(X2±X3)5*"^(F2±51s) 

and noting that the upper bound for EPR correlations is then 4. We see that this system 
demonstrates entanglement and the EPR paradox for a wide range of pumping strengths, 
these correlations persisting well into the region where the output fields are relatively 
intense and truly macroscopic. 

6. Entanglement via three concurrent nonlinearities 

We now turn our attention to a process in which the entanglement is produced in 
a single nonlinear interaction which combines three concurrent nonlinearities. The 
Hamiltonian we will investigate in this section is derived from work by Pfister et al 
[22], who raised the possibility of concurrent parametric down conversion in a single 
optical parametric amplifier (OPA). They also gave some solutions for equations of 
motion derived directly from the interaction Hamiltonian in the undepleted pump 
approximation, as well as experimentally observing triply coincident nonlinearities in 
periodically poled KTiOP04 

6.1. Properties of the Hamiltonian 

In this section we will also numerically solve the full equations of motion derived from 
the interaction Hamiltonian, as it is not possible to solve these analytically and it 
is known that the approximate analytic solutions and the full quantum solutions do 
not agree for arbitrary interaction strength in x^^^ systems |21 123 I2n]- We stress 
here that this analysis is not designed to give a full description of travelling-wave 
optical parametric amplification, which is not adequately described by our formalism 
(see, for example, Raymer et al 122!). However, this approach does let us examine 
the entanglement properties of the interaction Hamiltonian which we will later use to 
describe the interactions inside an optical cavity. A similar approach has been used 
previously to analyse, for example, the limits to squeezing and phase information in the 
parametric amplifier [2H1- 

A schematic of the nonlinear interaction is given in figure El showing the three 
inputs, which interact with the crystal to produce three output beams at frequencies 
tuo, oji and which may be equal. The interactions are selected to couple distinct 
polarisations, and the scheme relies on tuning the field strengths in order to compensate 
for differences in the susceptibilities since it is usually the case that Xyzy 7^ Xzzz- Note 
that X is the axis of propagation within the crystal. The mode described by hi is pumped 
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Figure 5. Schematic of a system described by the concurrent triple nonlinearity 
interaction Hamiltonian. 



at frequency and polarisation (u;o + uji,y) to produce the modes described by Si {ojq, z) 
and 0,2 {i^i, y), the mode described by 62 is pumped at {lji + lj2, y) to produce the modes 
described by 02 and 03 {uj2, z), while the mode described by 63 is pumped at (2u;i, z) to 
produce the modes described by ai and 02- 

The interaction Hamiltonian for the six-mode system is then 

ifint = {xSia\a\ + X2^24«3 + X3^3aI4) + h-c, (53) 

with the Xj representing the effective nonlinearities. In what follows, we will set Xj = X 
and the high frequency input intensities as equal, with vacuum inputs at the lower 
frequencies. These are the conditions that we have previously found to give maximum 
violation of the entanglement inequalities. 

Our first step is to analytically calculate the appropriate correlations in the 
undepleted pumps approximation. Setting ^ = Xi{bi{0)) as a real parameter, we may 
solve the resulting linear Heisenberg equations to find 



ai(t) = 


1 
3 


Aai(O) 


+ Bal{0) 


+ Ca2(0) 


+ ^4(0) 


+ Cds{o) 


+ d4(o) 


alit) = 


1 
3 


5ai(0) 


+ Ad[{0) 


+ 1^02(0) 


+ c4(o) 


+ DdsiO) 


+ c4(o) 


a2{t) = 


1 
3 


Cai(O) 


+ Dd\{0) 


+ ^02(0) 


+ Bdl{0) 


+ Ca3(0) 


+ ^4(0) 


al{t) = 


1 
3 


Dai{0) 


+ Cd{{o) 


+ Bd2{0) 


+ ^4(0) 


+ Dd3{0) 


+ c4(o) 


a^{t) = 


1 
3 


Cai(O) 


+ Ddl{0) 


+ Ca2(0) 


+ d4(o) 


+ ^03(0) 


+ 54(0) 


alit) = 


1 
3 


Dai{0) 


+ Cd\{o) 


+ Dd2{0) 


+ <^4(o) 


+ ^03(0) 


+ a4(o) 



where 

A = cosh2^t + 2cosh^t, 
B = sinh2^t - 2sinh^t, 
C = cosh 2^t — cosh ^t, 

D = sinh ^t + sinh 2^t. (55) 
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These solutions allow us to find the quadrature variances and covariances, 
V{K) = \[{A + Bf + 2{C + Df], 

V{Y^,) = \[{A-Bf + 2{C-Df], 

V{Xa,,Xa,) = ^-{C + D){C + D + 2A + 2B) , 

V{Y,,X,) = \{C-D){C -D + 2A-2B). 

For the van Loock-Furusawa correlations, written in shorthand as V3 since all three are 
equal, this then gives us 

1^3 = 1 [5(^2 + 5^) + %B{D - 2C) + 2A{AC ^ 8D ~ B) + U{C^ + D^) - 20CD] . (56) 
9 

The analytic expressions for the two types of three-mode EPR correlations (one- 
mode and two-mode inference) are rather more complex, therefore we will present the 
results graphically in figure El We note here that each correlation is identical for any 
permutations of the mode indices, due to the symmetries of the Hamiltonian. The 
two types of EPR correlations also have exactly the same shape, with the one shown, 
from (jSHj), being identical to the correlation of ()34|1 apart from a scaling factor of 4. 
Having seen that this system is potentially a good candidate, we will now integrate a 
more complete version of the equations of motion, taking into account depletion of the 
pumping fields. 

The full equations of motion resulting from this Hamiltonian are, we gain stress, not 
describing an optical parametric amplifier, but contain all the dynamical information 
which would be described by Heisenberg equations of motion derived from the 
Hamiltonian of . We now map the master equation derived from onto a Fokker- 
Planck equation for the positive-P function [21], making a correspondence between the 
operators dj, bi and the classical variables ai, We then find the appropriate stochastic 
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Figure 6. Approximate analytic solutions from the Heisenberg equations of motion 
for the van Loock-Furusawa correlations, V3, (any of the correlation functions of Q) 
and the EPR correlations of (|38|) . The value 4 represents the upper bound for true 
tripartite entanglement and a demonstration of the EPR paradox. 



differential equations in Ito calculus, 



dt 
<m 
dt 

<m 

dt 
dai 

~dt 



= -xa2a3, 
= -x«ia3, 



da2 
'dt 

~dt 

da. 

'di 

da 



dPt 

dt 

dPt 

dt 

dPt 
dt 



-XaJaJ, 



^ = X (/3i+a2 + f3tas) + r],{t) + r)e{t), 

- = X + Ptc^s) + ^5(t) + vm, 

r = X {P2at + P3at) + VW2 vm + VxFs 
^+ / / 
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Figure 7. Positivc-P solutions from (|57|l averaged over 11.92 x 10^ stochastic 
trajectories for the tripartite entanglement criteria and the three-mode EPR 
correlation. V3 is any of the correlation functions of Q and EPR is the correlation of 
()38|l . The line at 4 represents the upper bound for true tripartite entanglement. 



The complex Gaussian noise terms have the correlations 

r^=0, 7]*{t)7]k{t')=5,k5{t-t'). (57) 

An important point is that the pairs and {Pf) are complex conjugate only 

in the mean, due to the independence of the noise sources. This is necessary to allow 
the positive-P distribution to represent states which have a more singular distribution 
than the (5-function distribution of a coherent state in the normal P-distribution |3Uj . 
As always with the positive-P representation, averages of the variables converge to 
normally-ordered expectation values of the corresponding operators in the limit of a 
large number of stochastic trajectories. 

The above equations were numerically integrated and averaged over 11.92 x 10^ 
trajectories, with the results for the correlation functions of (j21) and (jHHj) being shown in 
figure The initial conditions were x = 10~^, Pj{0) = 10^ (real coherent states), and 
aj{0) = for j = 1,2,3. The horizontal axis is a scaled interaction time, ( = x|/5(0)|t. 
The minimum value of the correlations was approximately 0.02, which is very close to 
the zero predicted by the non-depleted theory. The van Loock-Furusawa inequalities 
are strongly violated over a range of interaction length, although phase noise from the 
depleted pumps eventually means that the violation vanishes, as it does with the EPR 
correlation. The EPR correlation for inferring one mode from the combination of the 
other two is again equal to that shown for one-mode inference, apart from the scaling 
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Figure 8. The spectral correlations for difTerent ratios of the pumping rate to the 
critical pumping rate, for the intracavity triply concurrent scheme. The parameter 
values are 7 = 10, k = 1 and x = 10~^. 

factor of 4. For the parameters we have used, everything is symmetric, with correlations 
not changing under changes of the mode indices. 

6.2. Intracavity 

We will now examine the experimentally more realistic case where the interaction takes 
place inside a pumped Fabry-Perot cavity, as previously investigated by Bradley et al 
IHj. In the best case, where the three pumping inputs and nonlinearities are equal, 
relatively simple analytic expressions can be found for the output spectral correlations 
equivalent to Q by following the usual linearised fluctuation analysis procedure. We 
make the proviso that these are not valid in the immediate region of the oscillation 
threshold, which occurs at a pump amplitude of e*'^ = 7k/2x, where 7 is the cavity 
damping rate at the high frequencies and k is the low-frequency damping rate. Below 
threshold, the steady-state solutions for the aj are all zero, while /^l** = e/j. Above 
threshold these solutions become = k/2x and aj'^ = ^/ {e — e^^)/x- We note here 
that, due to the presence of the square-root, there is an ambiguity in the sign of these 
solutions. However, closer analysis shows that all must have the same sign, whether 
this is positive or negative. The full spectral correlations have been presented in jS] for 
the damping ratio 7/K = 10 in which case the maximum violation of the van Loock- 
Furusawa inequalities is found at zero frequency. These zero-frequency correlations are 
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Figure 9. The spectral EPR correlations for different ratios of the pumping rate to 
the critical pumping rate, for the intracavity triply concurrent scheme. The parameter 
values are 7 = 10, k = 1 and x = 10"^. 



then found as 



4x2e2 (k7 + xe) 



When the damping ratio is changed so that K/7 = 10, the below threshhold resuhs are 
unchanged, but above threshold the spectra bifurcate so that the maximum violation 
is found at non-zero frequencies, as shown in figure ^| In figure ^1 we present the 
minimum of the correlation functions (maximum violation of the inequalities) for each 
of these cases. It is immediately obvious that the entanglement persists much further 
above threshhold in the bifurcated case. As this is no longer close to zero frequency, 
where technical noise can be a real problem, this may be a real operational advantage. 

We note here that these solutions give a limiting value of 2/9 at threshold and that, 
even though it does not result from a valid analysis, this threshold value serves as an 
absolute minimum which can be expected for these correlations, in the sense that 

4 > %(0) > 2/9 (59) 

represents the region where tripartite entanglement is shown for this system. We also 
note here that, above threshold, the output modes are macroscopically occupied, with 
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Figure 10. The spectral correlations for different ratios of the pumping rate to the 
critical pumping rate, for the intracavity triply concurrent scheme. The parameter 
values are 7 = 1, k = 10 and x = 10~^. 

intensities = (e — e**^)/x, so that, especially for the case with k > 7, genuine 

continuous variable tripartite is potentially available with intense outputs. 

7. Conclusions 

We have examined two different interaction schemes in terms of their potential as sources 
of continuous variable tripartite entanglement, in terms of both the well-known van 
Loock-Furusawa correlations and two three-mode EPR criteria which we have developed. 
While both give broadly similar results, the EPR criteria, which are written in product 
form, may be formulated more generally than the van Loock-Furusawa criteria which 
generally depend on knowing the correct combinations of the quadratures involved. In 
the specific cases we have examined there is a symmetry which makes a simple choice 
possible, but this will not always be the case. As for the actual schemes, we have shown 
that the one which mixes the outputs of three OPOs on beamsplitters and that in 
which the three entangled modes are created in the one intracavity nonlinear material 
have similar performance except in the far above threshold region. Which of these 
two schemes is preferable for practical purposes would seem to depend more on the 
robustness of the experimental setup rather than any inherent advantages that either 
may have. On the one hand, individual OPOs are familiar technology while the type 
of crystal needed for the concurrent scheme is relatively new technology. On the other 
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Figure 11. The spectral EPR correlations correlations for different ratios of the 
pumping rate to the critical pumping rate, for the intracavity triply concurrent scheme. 
The parameter values are 7 = 1, fc = 10 and x = 10~^. 



hand, it may prove easier to stabilise one cavity rather than having to simultaneously 
stabilise and synchronise three OPOs. 
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